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Kinetic mixing between the metric and scalar degrees of freedom is an essential ingredient in 
contemporary scalar-tensor theories. This often makes hard to understand their physical content, 
especially when derivative mixing is present, as it is the case for Horndeski action. In this work we 
develop a method that allows to write a Ricci curvature-free scalar field equation and discuss some 
of the advantages of such rephrasing in the study of stability issues in the presence of matter, the 
existence of an Einstein frame and the generalization of the disformal screening mechanism. For 
quartic Horndeski theories, such procedure leaves, in general, a residual coupling to curvature, given 
by the Weyl tensor. This gives rise to a binary classification of scalar-tensor theories into stirred 
theories, for which the curvature can be substituted for, and shaken theories for which a residual 
coupling to curvature remains. Quite remarkably, we have found that generalized DBI Galileons 
belong to the first class. Finally, we discuss kinetic mixing in quintic theories for which non-linear 
mixing terms appears and in the recently proposed theories beyond Horndeski which display a novel 
form of kinetic mixing, in which the field equation is sourced by derivatives of the energy-momentum 
tensor. 

PACS numbers: 04.50.Kd, 98.80.Cq, 95.36.-fx, 98.80.-k 


I. INTRODUCTION 

The evidence for accelerated stages of expansion in our universe’s cosmological evolution has been increasing over 
the last years. The first of such periods, cosmic inflation, would have occurred in the very early universe, giving 
birth to a flat, homogeneous and isotropic space-time, filled with thermalized radiation and imprinted with nearly 
scale-invariant and adiabatic perturbations. The second, late-time cosmic acceleration, reflects a much lower energy 
scale and only unfolds in the low redshift universe. Support for these phases of cosmological evolution comes from 
a number of complementary and increasingly precise probes that explore both the late and the early universe IDEO- 
The implications of cosmic acceleration for fundamental physics will be further scrutinized with the next generation 
of experiments, like the Euclid satellite, the Dark Energy Spectroscopic Instrument and the Square Kilometer Array 

PE). 

The simplest known mechanism for acceleration, a cosmological constant, can not satisfactory explain cosmic 
inflation without the introduction of a mechanism for its decay and a departure from scale invariance of the initial 
perturbations. Since inflation has to end, a more natural explanation is that it is produced by an additional dynamical 
degree of freedom, whose energy density eventually decays into dark matter and standard model particles Pi- For 
what concerns the late-time acceleration, the cosmological constant is able to explain current observations for cosmic 
acceleration, but it suffers of many theoretical problems, which make the investigation of alternatives a compelling task 
[9HTTj. A simple possibility to address the two phases of cosmic acceleration is the introduction of scalar degree(s) of 
freedom. Indeed, scalar fields are compatible with the symmetries of the cosmological space-time, can easily produce 
cosmic acceleration and occur naturally as limits of high energy theories of gravity. The search for models able to 
explain cosmic acceleration has triggered considerable interest in alternative gravitational theories (See [T2j for a 
recent review). 

In particular, scalar-tensor (ST) theories of gravity have existed in the literature since the early sixties, when 
alternative theories were developed in parallel to increasingly precise tests of gravity in the Solar System IT5] . The in¬ 
teractions present in old-school Jordan-Brans-Dicke (JBD) theories 114j constitute the first generation of ST theories 
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which was soon developed to a consistent framework for alternatives to Einstein gravity [15 , [H]- Recent develop¬ 
ments in extra dimensions and massive gravity have also uncovered new theoretical frameworks which produce viable 
modifications of gravity: the generalization of the interactions found in the 5-dimensional Dvali-Gabadadze-Porrati 
braneworld model HU led to the proposal of Galileon ST field theories m , which also arise naturally in the recently 
proposed dRGT ghost-free massive gravity m and bigravity [50] in the limit in which gravity decouples (see [5T i 55 
for reviews). While the aforementioned theories describe essentially different infrared physics, they are characterized 
by the same set of derivative interactions of the scalar field in the decoupling limit. 

The generalization of these interactions to curved space-time [53] naturally leads to Horndeski’s theory [53], which 
was first proposed in the early seventies. This is the most general ST action in four dimensions whose variation 
produces second order equations of motion. It characterizes the second generation of ST theories and encompasses a 
large set of models proposed over the past years. As such, this theory has attracted considerable attention as a way to 
unify ST theories and study their phenomenology as applied to late time cosmology [55H57], inflation [55], and local 
gravity tests [29] , Although this is the most general action whose variation produces equations of motion which are 
at most second order in derivatives, the completeness of Horndeski’s theory as the master framework encompassing 
all viable ST theories has been recently challenged: examples of theories beyond Horndeski indicate the existence 
of a third generation of healthy ST theories [50H32) . Given this large number of different models and the need to 
check their viability, many attempts have been made in formulating model independent observables as to test general 
properties of ST theories [25] [57] (33], as well as in the search for connections between apparently different theories 

m- 


This work provides a general procedure to investigate the properties of Horndeski’s theory by using, as a starting 
point, the way the scalar and the metric degrees of freedom interact. In fact, the aforementioned theories exhibit 
different degrees of kinetic interaction between the scalar and tensor degrees of freedom, a phenomenon known as 
kinetic mixing or kinetic braiding [35] . This property entangles the derivatives of the scalar and tensor held in the 
equations of motion, in a way that is unique to each theory. As we move from the simpler, JBD theories further into 
theories belonging to the second and third generation, the kinetic mixing becomes more intricate and new coupling 
structures appear. The idea is then to use the metric equations to remove all the instances of the curvature (which 
contains the second derivatives of the metric held) from the scalar held equation of motion. This procedure for 
covariant debraiding greatly simplihes the study of the properties of the scalar degree of freedom. 

We first illustrate how this procedure works for JBD ST theories and for cubic Galileons 135] and subsequently we 
introduce the debraiding procedure for general quartic Horndeski theories. This case represents a completely new 
situation, since couplings between curvature and scalar derivatives appear at the level of the action, and we found 
several new aspects for this extension that are not present in simpler cases. In particular, this procedure can not be 
completed in general due to a coupling between the scalar held and the Weyl curvature tensor (the traceless part of the 
Riemann tensor), which is not algebraically determined by the metric equations of motion. The debraiding procedure 
can also introduce spurious solutions to the equations of motion. However, we show that it is always possible to 
select the physical branch of solutions. We also found that in a specihc subset of quartic theories, Dirac-Born-Infeld 
(DBI)-like theories, the debraiding procedure can be performed in an exact manner, and both the Weyl tensor and the 
spurious solutions are automatically eliminated. This result is expected, as such theories are equivalent to Einstein 
gravity via a disformal redefinition of the metric. These two different behaviours under the debraiding procedure 
suggest to classify ST theories according to the possibility to remove all curvature ( stirred ) or not (shaken). 

The unbraided equations provide a new look into the properties of Horndeski theories: It unambiguously shows 
the couplings of the scalar field and sheds light to its behaviour within matter, as schematically shown in table [Tj 
Using our formalism we show that the DBI-like theories can present gradient instabilities in a radiation dominated 
universe with a sufficiently high pressure density, posing a serious challenge for the simplest among such theories. 
This instability can be easily avoided in non-DBI-like theories, for which the richer mixing structure can prevent the 
gradient instability. Finally, the unbraided equations allow one to generalize the disformal screening mechanism for 
scalar modifications of gravity (3SU323- This effect, which has been studied only for DBI-like theories in the Einstein 
frame allows the scalar field to evolve independently of the energy density if the scalar’s time evolution is non-negligible 
and the energy density of matter is sufficiently large. In this work we show that this mechanism is present in a broader 
class of quartic Horndeski theories, in which it may work under a more relaxed set of assumptions. 

We also show how the mixing structure acquires even more involved forms beyond quartic Horndeski theories. In 
quintic Horndeski theories the field equation has terms which are non-linear in the curvature, in the form of the Gauss- 
Bonnet scalar. Such non-linear mixing terms contain the square of the Weyl tensor and can not be debraided using 
our techniques. Theories beyond Horndeski introduce an even more subtle form of kinetic mixing, in which the scalar 
held is sourced by derivatives of the energy-momentum tensor (which can not be replaced using energy-momentum 
conservation in a covariant way). 

The paper is organized as follows. In section[TT]we will review the concept of kinetic mixing and covariant debraiding 
in old school and cubic ST theories. This study is extended to quartic theories in section III where we identify the 
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Old-School 

Horndeski 

3 rd generation 

Examples 

JBD /f(R) 

covariant Galileons 

covariantized Galileons 

Ostrogradski Degeneracy No d 2 (f> in £ 

2 nd -order eqs. 

Implicit constraints 

Kinetic Mixing 

algebraic 

T 

algebraic 

derivative 

V(T^'> 1 '), vt,... 

Formal Invariance 

under metric redefinitions 


C{(j>)g^ + D ((/>)(/>,, 

D(X, H 


TABLE I: Three generations of ST theories and some of their theoretical properties. The first line describes how the Ostrogradski 
degeneracy [38] is avoided as we move towards more complex ST theories. In the second line, the typical scalar-matter couplings 
that appear after the debraiding procedure is carried out. Finally the third line, shows which metric redefinition leaves the 
action formally invariant, amounting to a redefinition of the functions that specify the model. 


mixing structures, apply the covariant debraiding program and comment on the new subtleties that appear. In section 
|IV| we will explore the consequences of such mixing, both in the general case and specializing to the specific case of 
quartic DBI Galileon [23]. In particular, we will discuss the relation between the coupling to the Weyl tensor and the 
existence of an Einstein frame, the stability of the theory and the generalization of the disformal screening mechanism. 
In section [V] we will extend some of the discussion to quintic theories and to theories beyond Horndeski, showing how 
new matter-scalar couplings appears in such models. Finally, in section |VT| we will discuss our results and draw the 
conclusions. 

We will work in four space-time dimensions, use a —\- ++ convention for the metric, set the speed of light and the 
reduced Planck constant to unity h = c = 1. Summation over repeated indices is assumed. 


II. KINETIC MIXING: DEFINITION AND SIMPLE CASES 

In the introduction we have pointed out how contemporary ST theories can show a very complex mixing of their 
degrees of freedom in a way that is specific to the theory at hand. In particular, we have stressed how the couplings 
between metric derivatives and scalar field in the action will lead to second derivatives of one field to act as source 
for the other and to couplings between them, resulting in a very complex coupled dynamical system. In general, 
this mixing will not only complicate the numerical solution of the equations, but will also obscure their physical 
interpretation. 

One interesting way to simplify the scalar field equations is what we will call the covariant debrading procedure. This 
basically amounts to identify the scalar-curvature couplings that appear in the scalar equation, and use contractions 
of the metric equations with the scalar field derivatives to trade those for terms which depend on the scalar and matter 
fields. The outcome of this procedure will be a scalar equation of motion that depends on the scalar and matter fields 
and whose only second derivatives are those of the scalar field. 

There are several reasons for pursuing this idea. First of all, an unmixed equation of motion for the scalar field makes 
clear the interaction structure of the theory and the couplings between matter and the scalar field. This allows to use 
the debraided equation to study the stability of a given theory without the need to take into account also the metric 
equations. The kinetic mixing properties of a theory also determine its phenomenology; for example, kinetic mixing 
is necessary for any model to have a variable effective gravitational constant in cosmological scenarios BZl- Another 
reason is that this procedure allows a classification of different models depending on their matter and self-interactions. 
In fact, we will see that these are unique features of any given model and can help to distinguish between theories 
that are not equivalent via field redefinitions. Finally, we notice that this procedure has the additional advantage of 
being fully non-linear and covariant (not relying on a specific expansion or choice of background) and of using the 
Jordan frame matter stress-energy tensor which is covariantly conserved. This is simpler than rewriting the theory in 
the Einstein frame, which is not possible in general and leads to the energy-momentum being sourced by the scalar 
explicitly. 

In this paper we will mainly work with Horndeski’s theory [241 in its modern formulation [39] , described by the 









4 


following action 


Sn\g^, 4>] 


1 


d A xV^Y, 

i =2 


r 


(i) 


where 


A = G 2 (X,</>), (2) 

A = G s (X,</>)n<l>, (3) 

A = G 4 (x, cfyR + G 4 ,.y ((n^) 2 - (</>;^) 2 ) , (4) 

A = G 5 (X, $)G^ V - ^G 5 , x ((D^) 3 - 3 + 2(AA 3 ) , (5) 

are the quadratic, cubic, quaHic and quintic Lagrangian respectively. 1 Here X = — and Ocj) = , 

{<!>■, P v) n = ^ ; “”ai''' f denote contractions of the field’s second derivatives. As we have discussed in the 

introduction this theory has attracted considerable attention in recent years as a way to unify ST theories and study 
their phenomenology and hence represents the best framework for investigating the debrading procedure. 

We will introduce the essential features of kinetic mixing by presenting results for JBD theories and cubic Galileons 
(the simplest of Horndeski’s theories) in which all the basic features are already present. Quartic Horndeski theories 
will be presented separately in section m and |IV| while the novel kinetic mixing features introduced in quintic and 
non-Horndeski theories will be briefly discussed in section |V| 


A. Old-school scalar-tensor theories 


Let us start exploring the issue of kinetic mixing in JBD theories of gravity. Here we focus on a theory described 
by a coupling between the field and the Ricci scalar, a canonical kinetic term and a potential for the field: 2 

Sjbd = J d 4 x^g (^C(0) 2 R + X- V(<j>) + C m j , (6) 

The dynamics of the above theory is described by the metric equation 

C 2 G „„ + (g^nC 2 - Cl„) = ± ( T + g pv (X - V)) , (7) 


where the energy-momentum tensor is defined as T pu = ~^= 


8gP v ’ 


and the scalar field equation 


□ 0 - V' + CC’M 2 R = 0. (8) 

The kinetic mixing is reflected in the fact that the kinetic terms of the metric (~ R ^„) and of the held ~ VVC) 
appear in both equations. It is possible to debraid the scalar held equation by taking the trace of the metric equations 
and substituting it in Q 

The debraided held equation has the following structure 


(1 + m 2 Pl C’ 2 ) + 4^V - 2X^ (6 M 2 Pl {C’ 2 + CC A 


c. 


c 


+ 1 )= -£t 


renormalized kinetic term 


additional terms 


matter coupling 


(9) 


and contains no second derivatives of the metric. This simple example already reveals some of the debraiding features 
that will occur on more general theories: 


1 These names have historic origin. They refer to the power of the field in flat-space Galileons, for which Gi oc </>, ^ 2,^3 oc X , 
G 4 , G 5 oc X 2 . 

2 This is not the most general formulation of an old-school ST theory but such theory can always be mapped to this form with suitable 

field redefinitions Bsmn 







5 


1. There is an explicit coupling to matter, leading to an environment-dependent effective potential. In this simple 
case it is proportional to the trace of the energy-momentum tensor T, as could be anticipated from coupling 
between the Ricci scalar and the scalar field function G in eq. ([ 8 ]). Notice that this coupling is one way , in the 
sense that a minimally coupled matter source will still have the matter stress-energy tensor conserved. 

2. The coefficient of the second derivative term gets renormalized by a function of the field, showing how kinetic 
mixing can affect the stability properties of the scalar field equation. In this case the coefficient is strictly 
positive and hence no instabilities can be dynamically generated. However, this will not generally be true for 
more complex theories. 

3. New terms not involving second derivatives appear in the equation, coming from the contraction of the first 
derivative terms in the metric equations. In particular, the potential term for the field is modified and a new 
term involving first derivatives of the field appears. 

It is worth stressing that matter is minimally coupled and therefore the energy-momentum tensor is covariantly 
conserved = 0. This would not be true if similar results were obtained by expressing the theory in the 

Einstein frame by a redefinition of the metric. 


B. Cubic Horndeski theories 


Cubic Horndeski theories ^ are characterized G 4 = M‘pj2, G 5 = 0 and generic functions G 3 and G 2 . As we will 
see, they contain richer forms of kinetic mixing in curved space-time. Here we will explore their behaviour for the 
simplest non-trivial example, the cubic Galileon: 


Scg = 


d 4 


M 2 


X 


X V ~9 ( + X + -^3 □</> 


( 10 ) 


The second derivatives of the field present in the last term produce a coupling with the affine connection, which in 
turn introduces a term involving the curvature in the field equation 


s s 

— = A^ 3 </> ^R^tj) v + terms without curvature, 
dtp 


( 11 ) 


(one can alternatively see the the emergence of the Ricci tensor through the anti-commutation of covariant derivatives, 
which appear anti-symmetrically in the equations of motion). Cubic Horndeski theories are known as kinetic gravity 
braiding (KGB) [35] [42] for this reason. 

Just as in the JBD case, it is possible to use contractions of the metric equations to solve for the curvature coupling 
in (111. The only difference is that one has to contract with both the metric and The resulting debraided field 

equation reads 




2X 2 . 
Mg A 6 ' 


4X 

Ml A 6 ' 




renormalized kinetic term 


4X 2 
Ml A 3 


extra terms 


1 

at - 


(□ 0) 2 - 


MIA 3 


+TX) = 0. 


( 12 ) 


higher derivative interactions 


coupling to matter 


We note the following features: 

1. The coupling to matter has two contributions, a conformal one (proportional to the trace T and weighted by 
X) and a disformal one given by the contraction of the energy-momentum tensor with tf>^tj> „. The disformal 
part is particularly interesting, as it indicates that radiation would have non-trivial effects on the field in this 
type of theories. 

2. The kinetic term is renormalized due to the braiding by both conformal oc g as well as disformal oc 
terms. Unlike in the old-school case, the corrections are not positive definite any more. 

3. There appears a Galileon term constructed out of anti-symmetric, non-linear second derivative contractions. 
This term is not renormalized in the unbraided form of the equations. This type of terms are responsible for 
the Vainshtein screening mechanism in cubic theories that allows these theories to fit local gravity tests. 
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The authors of (5HJ used the debraided equations to study causality and stability of KGB theories. In what follows we 
will extend the same program to quartic Horndeski theories and discuss the new features and subtleties that appear. 


III. KINETIC MIXING IN QUARTIC HORNDESKI THEORIES 


In this section we extend the debraiding formalism to Quartic Horndeski theories defined by the fixing Gs((/>, X) = 0 
and Gs((p,X) = 0 while leaving the other two functions arbitrary. 3 

These theories, defined in Q, introduce a range of new kinetic mixing terms which arise in the field equation via 
both the cancellation of higher derivatives (due to the anti-symmetric structure, as in cubic theories) and the direct 
coupling between the field and the Ricci scalar in the action: G 4 (X, <f>)R. 4 The metric and scalar field equations are 
reported in appendix [Aj while here we focus on the mixing terms which appear in the equation of motion with the 
following structure: 

^ = -2G 4 , A G“^ iQj8 + 2G 4 ,* a -0’> /3 (^2 <f>. a ’ x R X p + - n <^/3 - f ^ 

+G 4i ^.R — 2G^x<t>Raf3<t>’ a <t>’^ + terms without curvature . (13) 


We note two distinguishing features that did not appear in lower order theories. First, there are second derivatives 
of the scalar field multiplied by Ricci curvature. Therefore, using the metric equations will introduce new products 
of field derivatives into the equations, which will in general fail to be linear in the second time derivatives and hence 
may introduce spurious solutions to the equations of motion. Second, the derivatives of the field also couple to the 
full Riemann tensor. This term can be rewritten in terms of the Weyl tensor as 

P p ua0 ^Vfj.vafU T 9 p [aP/3]i/ ^P*9 P [a9 0]v 3 (14) 


where the square brackets stands for antisymmetrization of the n encompassed indices with weight 1/ro! . This is a 
necessary step in order to split the Riemann tensor into its trace part, solvable from the metric equations and into 
its traceless part, which cannot be solved for using contractions of the metric equations. Besides introducing new 
interesting features, both aspects represent an obstruction to the debraiding process. The spurious solutions are a 
technical complication that can be surpassed, as we will explain in section [Til A| On the other hand, the Weyl coupling 
is not devoid of physical meaning and its consequences are explored in section [IV A 

In the case under investigation, where G 4 = G 4 (</>, X) and G 2 = G2(0, X), the field Euler-Lagrange equation can 
be written in the following compact form 

L pv (j>^ v + V + + Q pvapp °^^ a 0^ P a = 0 , (15) 


with 


D iV 


p P va0 


XX RJ 

2 (G 4 , a + G 4 , aa ) G pv + 2G 4 , AA (cj>’ a G a p ^ a v p ) , 

(3G 4 , 0A - 2 G 4 , 0AA ) {g pu g a P - 9 pa g vf) ) + 2G 4 ,^ AA (2^0>- 2 ^rg a ^) , 


G 2 , a + g 4 , aa (g^^ 


-XR 
3 


— I G 2 , AA — 0^4, 


QlAlSOifipCT 


V 


G 4 , aa {g^g a0 g pa - 3 g pv g ap g^ + 2^VV p ) 

G 4 , aaa (2 ^g^g^r - 2 ^’W^ ,/3 + W - 9 pa g a0 )) , 

G2,(j> — 2 XG 2 ^ x — 4G 4i ^ A (/>’ q G q;/ 3(()’ /3 + Gi^R. 


(16) 

(17) 

(18) 
(19) 


where the first term is linear in second derivatives, the second one is a potential term that depends at most on first 
derivatives of the field while the last two terms, despite being respectively quadratic and cubic in derivatives, are 
linear in second time derivatives and non-linear only in mixed spatial derivatives. 


3 While the second choice is just for practical use, as it simplifies the calculations without significantly altering the analysis, the first one 
is more important. In fact, taking Gs((f), X) 7 ^ 0 would introduce significant deviations as we will discuss later on. 

4 Note that the dependence of X generates derivatives of R in the field equation, which cancel exactly with counter terms stemming 
from the pure field part m ■ This follows from the defining property of Horndeski theories, namely the absence of higher than second 
derivatives in the Euler-Lagrange variation. 
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As it has been done in the previous cases we can eliminate the curvature-field couplings with suitable contractions 
of the metric equations with scalar field derivatives. In this case the structure is more involved and hence we describe 
it in a schematic way, leaving the full expressions to appendix [A] The debraided field equation is 

+ V+Q t T+ Q {T )t' 0 T a ^ a + + K pva ^ 

+ + Q pvapp °) = o, (20) 

with 


L pv = ( g 0 + Q(T)^T a ^ a + Q t T ) g pv + (So + S {T )^T a ^ a + S t T) 

+ C {t) ^T^^ +C t T» v + C w <t>' a ^W» a v p, ( 21 ) 

P^ aP = v 4s (g^g af) - g pa g vf) ) + V 4D (</>’'W /3 - ^’V^) , (22) 

K>* vaP = W D 2 + 4Al V"s a/3 + 4X ^^ v g a P) , (23) 

Q^ pa = V 5g (g pv g afi g pa - 3 g pv g ap g^ + 2 g»°g va <f p ) 

- V 5 .Y (2 rg va g &p r - 2 <S>' p g va g p °^ + W (g p °g af) - g pa g^)) , (24) 

H pv a [tp* = Wi ((g^g°<p -^V' 9 ) (p P p a +2Xg Pa ) 

+ 3 (<K0’V“ - grw* - X (g pv g afi - g pa g v0 ))) g p ” 

+ W 2 - g pu 4>’ a (j)’^ - X (g^g^ - g pa g vf) )) (^^ + 2Xg pa ) . (25) 


Here the coefficients V,Gi,Si,Ci, V;, Wj, <2* depend on </> and X through G 4 , G 2 and their partial derivatives and are 
fixed once a specific model is chosen. In table In] we schematically report the coefficient structure for three models of 


quartic Horndeski, while the general expressions can be seen in appendix A 2 


After the debraiding process some of the terms in the debraided equation take a form analogous to the unbraided one 
but with “renormalized” structure coefficients. Extra terms which couple matter to the scalar field are also introduced. 
In this regard, it is interesting to note that these coefficients (see appendix [A]) have a common denominator structure 


~ (G 4 - 2IG 4 ,x) _ " (G 4 - IG 4 ,x)- m . (26) 

In particular, the coefficients can become singular for certain values of the field. However, the first factor is inversely 
proportional to the effective Planck mass, defined as the coefficient of the second time derivative in the graviton 
propagation equation. In homogeneous and isotropic backgrounds, and restricting to quartic theories, this has been 
shown to be m 


Ml ff = 2(G 4 - 2XG 4 ,. y ) . 


(27) 


This equation can be covariantized right away, suggesting its validity on general backgrounds. This dependence is 
consistent with the fact that the debraiding procedure substitutes the curvature terms with the energy-momentum 
tensor via the metric equations and hence it suppresses the new terms by a Planck mass factor. Therefore the 
singularity of the unbraided equations is related to a physical singularity, as the coefficient of the graviton kinetic 
term vanishes. Moreover, requiring the positivity of the effective Plank mass also implies that the other factor in the 
denominator will never become singular as long as G 4 is positive. 

The first term in equation (20), as can be seen from (21), is linear in second derivatives and is of particular interest 


as it contains the couplings between matter fields and second derivatives of the scalar field. This represents a new 
type of contributions with respect to previous cases and is important for the stability of these models. As we will 
discuss below, depending on the nature and evolution of matter fields, instabilities may occur for certain theories in 
the presence of matter. We point out again the presence of the Weyl term, indicating that there is still a residual 
dependence on the second derivatives of the metric whose physical meaning will be discussed in section IV A 


The Qi terms contains other couplings between the scalar held and matter which involve at most first derivatives 
of the scalar held, while the other terms are non-linear derivative interactions, responsible, for example, for the 
Vainishtein screening mechanism. However, among those, the terms proportional to H i can be quadratic 

in second time derivatives on non-trivial backgrounds. This dependence may introduce spurious solutions, which can 
nonetheless be distinguished from the physical ones (see section III A). 







-Gt(4>,X) 


Qt Q(t) St S(t) Ct C(t) Qt Q(t) Cw V4 d V4s V5 d V5 b W1 W2 Wd 2 


^1 - 2 AY A 4 (37l 


\/l ~ 2A{4>)X/X± 

(i+WA 4 )") 144] 

(1 + (X/A 4 ) 


(33 1 


0 0 
0 0 
/ / 


0 0/0 
0 0/0 
/ / / / 


0 

/ 

0 


0 

/ 

0 


0 0 0 0 

0 / / 0 

/ 0 0 / 


0 0 0 
0 0 0 
/ / / 


0 

0 

0 


TABLE II: Debraided coefficients scheme for three models of quartic Horndeski action. The ticks indicate which one of the 
debraided coefficient is present for each theory. 


A. Spurious solutions and their avoidance. 


The covariant debraiding procedure introduces unphysical solutions in the field equation. Contractions of the metric 
equations with 0 ;At „ and (necessary to the covariant debraiding process) lead to the introduction of quadratic 

powers of second time derivatives in eq. p0| . Therefore, two branches exist for the solutions, the physical one and a 
spurious mode, which has to be disregarded. 5 

The debraiding procedure is equivalent to summing the field equation and the combination 


M = ( c 0 g^ + c 2 ^ + c 3 <^ ;i/a </>,«) - T^) = 0, 


(28) 


where anc ^ c o~ c 3 are chosen to cancel the braiding terms in the field equation. The problematic 

c t) 2 terms appear in the structure term fjP va Pp a and in eq. (20) as a consequence of the contraction of 


second derivatives of the scalar field with the metric equations that also contain such terms that lacks of the required 
antisymmetric structure required to avoid non-linear terms. Note that the coefficient of 0 2 involves contractions of 
spatial derivatives and hence solutions on simple geometries (e.g. FRW) will not display the complications associated 
with the spurious solutions. 

We can use a trick to help us pick up the correct branch for the solution. Instead of adding A4 to the field equation, 
we can deal instead with 


SC 

— + eM = 0. 

dip 


(29) 


The above equation allows one to regard the debraiding procedure as a continuous deformation of the field equation 
depending on a parameter e: it is now possible to interpolate between the original field equation e = 0 (which is 
linear in second time derivatives) and the unbraided field equation e = 1. Once a time direction has been chosen (e.g. 
through an ADM decomposition or an explicit choice of coordinates), one can schematically study the field equation 


ecj) 2 + B(f) + C = 0 , 


(30) 


where e has been left explicit only in the term quadratic in second time derivatives. This form is guaranteed by the 
fact that all non-linearities in second time derivatives are introduced by the unbraiding terms, and are therefore linear 
in e. The existence of real solutions will be determined by the condition B 2 — 4 eC > 0. 6 In that case there will be 
two solutions 


4>=y £ (-B±VB^4eC) =^(l±l) T ^+0(e), 


(31) 


where the last equality entails an expansion on e. We can identify the unphysical branch as the one associated with 
the + sign, as it is not mapped to a solution of the original equation when e —> 0. As the physical solution should 
depend continuously on e. the solution associated to the minus sign is the physical one while the plus sign leads to a 
spurious solution which is not originaly present. 

Summarizing, one can choose the p hys ical solution of the unbraided equations after choosing a time coordinate by 
writing the field equation in the form (30), setting s = 1 and integrating in time (j> using the solution with the negative 
sign. 


5 A cartoon example of this: If our original equation has the form fix) = a, taking the square will introduce a spurious solution for x, 
corresponding to fix) = —a. We thank J. Beltran for helping us clarify the introduction of spurious solutions. 

6 It is theoretically possible that the debraiding equations have no real solution, while the original ones do (if B 2 < 4 C). In that case one 
would need to diagonalize the system of differential equations in some other way in order to numerically integrate the dynamics. 
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B. DBI-like Galileons 


In this section we are going to apply the debraiding method described in the previous section to a specific set of 
models for which the debraided equations are greatly simplified. We will then use these models as a reference when 
discussing the applications of the debrading method. Quite remarkably, these models are a generalization of the 
quartic Dirac-Born-Infeld (DBI) Galileon j33], in which 


Mi 


G^,X) = ^ ] Jl-2A^) A4 


X 


(32) 


where A is a new mass scale, while G 2 is left generic. 7 

The debraided field equation for this theory has the remarkably simple form: 


with 


+ V + QtT + Q {T )^T a ^’ a + = 0 


= g 0 g^ +So<P<l> v + C t T» v , 

p lu,af> = V 4B (g^g a0 - g^g^) + V 4D ~ ^9 aP ) , 


(33) 


(34) 

(35) 


where again the exact form of the coefficients can be found in the appendix A 2 a| Here we note some interesting 
features in the debraided field equation. First, the linear term contains only one coupling between and the scalar 
field second derivative, while two terms couple to the trace of the energy-momentum tensor and its contraction with 
first derivatives of the field (as an effective potential). Second, it only contains non-linear derivative interaction terms 
of order (VV<(>) 2 : the (VV<(>) 3 terms characteristic of quartic theories cancel in the debraiding procedure. Finally, 
both the Weyl coupling and the n on-linear second time derivatives cancel from the equations. Hence, the spurious 
solutions analysed in section III A are absent. 8 

The reason for this simplicity is that the theory can be casted in a much simpler form by means of a field redefinition, 
as the DBI Galileon is equivalent to Einstein gravity plus a disformal coupling to matter pa ns]. Because of this 
simplicity, DBI Galileons offer a toy example of kinetic mixing in more complicated Horndeski theories which also 
offers a viable and interesting alternatives to inflation m and as a mechanism for present day acceleration 

If we further restrict the class of models and consider a constant A in the definition of G 4 so that 


Mi I ~X 

G4 = ^Vl-2 F , (36) 

then the equation takes the even simpler form 

L^^ v + V = 0, (37) 

where 

^ = ( G ** + idSrAi)- ( G — - T*T^m) ^ 

I _I_ ^ TV”' 

A 4 1-2X/A 4 
V = G^^ — 2XG2,tj,x 

where only the coupling between matter and second derivatives of the field remains. Also notice that the non-linear 
derivative interaction terms have completelly vanished in this case. Thus this model does not have a Vainshtein like 
screening mechanism, althought it might possess a D-BIonic screening m■ This model is equivalent to the simplest 
disformally coupled theories. 


(38) 

(39) 


7 We could have added a field dependent Planck mass, but this would not modify most of the results. 

8 DBI theories have other interesting properties. For example, when regarded as effective theories they only require that the second and 
higher derivatives remain small <C A n+1 while they allow for arbitrary values of first derivatives [45j . 
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IV. CONSEQUENCES OF KINETIC MIXING 

In this section we explore some of the consequences of kinetic mixing, namely the coupling to the Weyl tensor, the 
stability in the pressence of matter and the possibility of screening modifications of gravity by kinetic mixing. The 
discussion will present both results for general quartic theories and their simpler DBI-like counterparts. 


A. Coupling to the Weyl tensor and existence of an Einstein frame 

One of the most salient features of quartic Horndeski theories, which is highlighted by our procedure, is the coupling 
between the Weyl tensor and the derivatives of the scalar field. The specific form of the coupling is 

> (40) 

where the right hand side identifies the electric part of the Weyl tensor with respect to the constant <j> hypersurfaces . 9 
Note that due to the symmetries of the Riemann tensor the above is the only independent contraction that one can 
form with derivatives of the scalar. By the traceless property of the Weyl tensor, it only couples to q L >- a p = 4>-, a p — H<t>g a p 
and X) iL/ = + 2Xg lil/ . 

The Weyl tensor is not determined algebraically from the metric equations. Instead, it obeys a propagation equation 


= V\m g »]P + g^V.jG 0 


(41) 


and it is therefore determined by the other fields through differential identities. One can generally decompose the 
curvature into a trace part (Ricci) and a traceless part (Weyl) using equation |l4| ). In Einstein’s theory, the Ricci 
curvature vanishes in the absence of matter. Therefore the Weyl tensor fully characterizes vacuum effects such as 
gravitational waves and tidal forces, both of them sourced by derivatives of the energy-momentum tensor (as given 


by the above propagation equation after substituting G^ —> 8ttGT^ u ). One can formally invert eq. (41) and write 


the Weyl tensor with a non-local dependence on the energy-momentum tensor. 

In quartic theories, however, the direct coupling to the Weyl tensor implies that the scalar field has a new interaction 
with space-time curvature. The effects of this coupling are more difficult to interpret, since in ST theories non-trivial 
configurations of the scalar field can also produce non-zero Ricci curvature in the absence of matter. This feature 
prevents the Weyl tensor from fully characterizing curvature in the absence of matter, making it harder to link the 
coupling (40) to vacuum effects such as tidal forces or gravity waves. However, in some situations the scalar field 
contribution to the metric equations might be subdominant with respect to matter, and we can recover the usual 
notion of the Weyl tensor describing such effects. At least in those cases, the Weyl coupling will affect the way in 
which the scalar field is sourced by the gravitational field produced by distant matter. 

Another important consequence of the occurrence of the Weyl tensor is that it obstructs the debraiding procedure, 
as it can not be obtained through contractions of the metric equations. This fact can be related to an important 
property of these kind of theories which is at the foundation of our classification into shaken and stirred theories. 
The presence of the Weyl tensor in the equations of motion can be related to the lack of existence of a (local) field 
redefinition which renders the kinetic term for gravity canonical, i.e. of the Einstein-Hilbert form. The diagram in 
figure |T] shows how field redefinitions in the action correspond to linear transformation of the equations of motion, 
where the transformation matrix is given by the Jacobian of the field redefinition m- The Weyl tensor represents 
an element that can not be “rotated away”, therefore indicating that the kinetic term for the metric (given by £ 4 ) 
can not be made canonical by a local field redefinition. Non-local redefinitions might be an exception, as they might 
“undo” the effects of the propagation equation for W flva p (41). 


Hence, the fact that DBI-like theories do not produce a coupling to the Weyl tensor is fully consistent with the 
existence of an Einstein frame for this class of theories. It is known that only a very special subclass of Horndeski 
theories can be related to Einstein gravity via a local field redefinition. In reference [35] this was investigated 
using general disformal transformation, showing that the most general quartic theory that can be mapped via a 
special disformal transformation to its Einstein frame version has to take the very special DBI-like form in which 
Gi((f>, X) = A((j>)y /1 — 2 B{cj))X. An argument as to why more general field redefinitions involving the scalar field will 
not accomplish this goal is given in appendix |B| 


9 This implicitly assumes that is time-like. We thank I. Sawicki for pointing out the electric character of the coupling. 
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FIG. 1: Field redefinitions and linear transformations of the equations of motion. 


B. Stability in the presence of matter 


An important application of the debraided field equation is to study the stability of the theory in the presence of 
matter. 


When considering DBI Galileons, equation (33) shows a flaw of the theory, as the kinetic mixing term might becomes 


problematic in the presence of large matter pressure. More precisely, if the energy-momentum tensor of matter contains 
a positive isotropic pressure term Xk„ D pg^, p > 0 , then the speed of sound of the field perturbations can become 


imaginary. For the simplest case ([36]) with G 2 = A, the evolution equation (371 reads 


(A 4 - X) n</> + = 0, 


(42) 


where 1 — 2 A/A 4 7 ^ 0 as it is related to the effective Planck mass (I27J) . One finds that the speed of sound squared 
(given by the coefficient of the second spatial derivatives) can become negative, leading to a gradient instability. The 
critical value of the pressure above which this happens is approximately 


p c ~ max(A 4 , A'), 


(43) 


where the estimate is obtained by looking at the sign of the ii component of eq. (42), where it has been assumed 
that A 4 > 0 in order to make the second time derivative coefficient positive when p = T ao is large in order to avoid 
ghosts. The occurence of such instability is equivalent to the failure of the coefficient of in (42) to have Lorentzian 
signature and hence of the field equation to be hyperbolic. 

DBI-like theories can therefore develop a gradient instability if the term proportional to the energy-momentum 
tensor dominates (see also l52| for some remarks in this direction). This problem is particularly acute if we are 
interested in theories able to explain cosmic acceleration, for which A needs to be a very low energy scale. In this case 
we may spoil the predictions of homogeneous cosmology at early times, particularly during radiation domination when 
p ~ p/3. Therefore, the simplest solution is to raise the value of A to make it higher than the reheating temperature . 10 
Another solution is to give A a field dependence (equivalent to the more general DBI-like theory (33)) so that the 
critical pressure is always larger than the cosmological one . 11 

The last possibility is the re-introduction of non-linear derivative self-interactions by modifying the cubic term such 
that G;i_x 7 ^ 0, without modifying G 4 . These terms would dominate when the spatial derivatives become large, and 
may act to stabilize the equation at some finite wave number, given by the condition 


p/A 4 ~fc 2 G 3 ,x<^. 


(44) 


This estimate relies on te scaling of the derivative self-interactions in a cubic theory, which are ~ G 3 ^(VV^cxfc 4 . 
Even if this modification stabilizes the perturbations, it would introduce large spatial gradients which might spoil the 
homogeneity and affect cosmological observables. 

Finally, it is possible to avoid the gradient instability by making a different choice of G 4 . As it was shown in section 


III theories different from the DBI Galileon have a richer mixing structure, as given by equation (15). The additional 


terms present in the general case may balance the ones leading to the gradient instability. For example, the general 


10 Raising A 4 beyond the sca le of reheating is not necessary: inflation requires a negative pressure and therefore leads to no gradient 
instability in equation j42[ . Ref. m considers an inflationary scenario in which A 4 is negative, but so large that the energy density 
never flips the sign of the kinetic term. 

11 The instability might be resolved by self-consistently accounting for the evolution of matter. However, even if the dynamics does not 
lead to singularities, it may cause large inhomogeneities incompatible with cosmological observations (see for an example). 
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debraided equations for a quartic theory contain a term for which the spatial and time derivatives have always 

the correct relative sign and no gradient instabilities occur. It is therefore possible to avoid the gradient instabilities 
by constructing a theory in which such terms are sufficiently large. Of course, in order to assess the viability of these 
models, a full dynamical analysis is required, but this goes beyond the scope of the present work and is left for further 
studies. 


C. Screening scalar forces by kinetic mixing 


Alternative theories of gravity typically introduce additional forces, which may alter the predictions in the local 
system and render them incompatible with local gravity tests (for a review, see ref. |54 | ). However, some theories 
provide screening mechanism , which hide the effects of scalar forces via non-linear interactions of the held. Although 
these mechanisms have been mostly studied by making these couplings explicit (e.g. working in the Einstein frame), 
they can also be identified in a minimally coupled description [55| . 

We conclude the discussion of the consequences of kinetic mixing by noting that the structure of the mixing terms 
allows to generalize the previously proposed disformal screening mechanism to a more general phenomenon, based on 
kinetic mixing properties and present in a larger class of theories. The disformal screening mechanism was introduced 
in the context of disformally coupled theories, which are the Einstein-frame version of (and therefore equivalent to) 
the DBI-like Galileons that we considered in section IIIB Its action is based on two observations, see eq (33) 


1. If the held is static (no time derivatives) and only disformally coupled (M p is constant) then the Held decouples 
from non-relativistic matter [56] . 

2. If the Held evolves in time and the energy density is non-relativistic (only T 00 = p contributes significantly), the 
field evolution becomes independent of the energy density (see refs. [36] [37] for details). This property does not 
rely on the specific form of the conformal and the disformal coupling. 

The efficiency of this mechanism to reconcile DBI-like theories with local gravity tests is difficult to investigate in 
practice, as it requires considering simultaneously spatial and time dependence (demanding that the time-evolution 
is a sub-dominant effect clearly spoils the existence of a screened solution 152. although numerical studies seem to 
confirm the screening effects [58]). Moreover, pure DBI-like theories have to face the issues related to stability in 
the presence of matter with non-negligible pressure, as described in section |IVB[ as well as other difficulties to fulfil 
laboratory tests [55] , 

Ultimately, the main requirement for the disformal screening mechanism is that the coefficients of depend on 
the energy-momentum tensor. The identification terms containing the matter energy-momentum tensor in the second 
derivatives of the scalar in the debraided field equation (16) indicates that the disformal screening mechanism, or 


variations thereof, can occur in a much larger class of quartic Horndeski theories. Among others, those terms include 
contributions proportional to 


TU<j>, ^T^Dcj,, Tcj) ,a 4>’P(j>. a p , (45) 

where only the first one is present in DBI-like theories and the disformal screening mechanism. Therefore, the richer 
braiding structure of quartic Horndeski theories has the potential to soften some of the requirements needed for 
the disformal screening and alleviate some of its problems. Because of its extended generalities beyond disformally 
coupled/DBI-like theories, we propose referring to this mechanism as screening by kinetic mixing. 


V. BEYOND QUARTIC THEORIES 

In this section we comment on the issue of kinetic mixing in more general ST theories. We will not go into the same 
level of details as for the quartic, limiting our analysis to the mixing terms and pointing out the novel structures that 
appear. 


A. Quintic Horndeski theories: Non-linear mixing 

In the case of quintic Horndeski theories we expect terms quadratic in the curvature to be present in the scalar 
field equation. Due to the second order nature of the theory we will only find terms that do not introduce higher 
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derivatives of degrees of freedom. One example of this is the Gauss-Bonnet term Q = R 2 — 4 R IMV R tlv + 
that would be generated in the field equation as 

£gb3/W^/'W, (46) 


since the above Lagrangian is equivalent to a particular choice of Horndeski functions with G§ t x 7^0 [28]. 

These theories hence produce a new form of non-linear mixing. Covariantly debraiding the field equation will 
therefore require (at the very least) contractions of the metric equations with the Riemman and the Ricci tensor, 
and then further contractions similar to those needed in section III Whether such terms can be debraided using 


the techniques introduced in the previous sections lies beyond the scope of this work, but at the very least we can 
anticipate the difficulties we already encountered for the quartic theories. In particular, the coupling to the Weyl 
tensor occurs already at the level of the action. This is not manifest from the simpler form shown in eq. (15]), but can 
be obtained in an equivalent form obtained by integrating by parts (e.g. in Horndeski’s original paper |zT| . where 
dual of the Riemann tensor appears in the action). 


B. Theories beyond Horndeski: derivative mixing 


The previous sections have shown how Horndeski theories feature a form of kinetic mixing that is algebraic in the 
energy-momentum-tensor. Healthy non-Horndeski theories HE EE] display a novel form of kinetic mixing, in which 
the energy-momentum tensor enters the scalar field equation through its derivatives. These theories have received 
attention recently, including studies in the context of late time acceleration [32], [60], inflation [61] and local gravity 
tests [62U64] , 

The simplest examples of theories beyond Horndeski are the ones originally proposed by Bekenstein [55]. These are 
ST theories with an Einstein- Hilbert term for gravity, a k-essence Lagrangian and a matter Lagrangian constructed 
out of a metric that explicitly involves the scalar field. 


Sb [g^uA,^] = / d 4 x 


' _M 2 _ __ _ N N 

v-9-^-R[9ap\ + v^gG 2 (x, (j>) + (<?/«/, VO 


The novel ingredient is that the matter Lagrangian is constructed using a general disformal metric 

9nv\g a pA] = o(x, 4>)9 I u' + D{x,4>)4> tlx ^> jV . 


(47) 


(48) 


When written in the Jordan frame via a non-trivial inversion of eq. (481, these theories have been shown to be 
non-equivalent to any Horndeski theory [30] [46] unless Cj[,D t x = 0 Indeed, their Euler-Lagrange variation yields 
equations with derivatives higher than second order. 

This seems to suggest that the theory propagates an extra degree of freedom. However, it was found that an implicit 
constraint exists in the equations of motion for the metric which allows to remove all the higher time derivatives from 
the equations of motion and cast the dynamical equations in a second order form. This procedure uses a contraction 
of the metric equations in a manner analogous to the procedure performed in section |III| and thus the formulation 
using implicit constraints (which is second order) introduces the energy-momentum tensor in the equations of motion. 
Here we sketch the basics, the interested reader is referred to ref. [30] and appendix [c] for further details. 12 

The novelty of this class of theories is that the terms introduced with this procedure involve derivatives of the 
energy-momentum tensor in the field equation, therefore providing a new form of kinetic mixing. This can be seen 
by considering Sb , <P- % lA■> the Jordan frame version of (47), in which the equations of motion are obtained for g a p 
(rather than g MJ/ ) using the inverse of the matter metric (48). The inverse relation between the metric has the same 
structure, g IJV [g a p, q i] = A(X, <p)g^ + B(X, where A, B can be obtained implicitly given the form of C, D. 

We note that it is possible (although cumbersome) to express the Jordan frame theory in terms of g a p. This is done 
for two simple cases in appendix [C] 

By using the chain rule and the properties of the Jacobian of the transformation between the metrics, it is possible 
to write the field equation without higher order field derivatives as 


(r K r 


- <5^V m + G^ (A^ 



= 0 , 


(49) 


12 See I66| for similar remarks in a more general setting and a generalization of Bekenstein’s program to find the most general ghost-free 
theories that are related to Einstein gravity by a field redefinition. Some of the generalizations proposed there have also been explored 
in ref. [67l . 
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where A,B define the inverse of (48) and all barred quantities are meant to be evaluated in terms of g^ v using this 
relation. The kinetic mixing factor is defined as 


/K Ml V g A-A 


,x 


X 


2B xX 2 


(50) 


where the energy-momentum tensor has the usual definition in terms of a variational derivative with respect to g^ u . 
This factor enters through the first term in the field equation (491 and thus introduces a derivative form of mixing 
between the matter and the held, even though the scalar and matter Helds are minimally coupled. Of course, it 
remains to be proved that the curvature stemming from the barred Einstein tensor G^ v can be traded for algebraic 
couplings to matter and the scalar Held in an analogous way as it has been done for the quartic Galileon. In appendix 
|C 1| it is shown how for the pure conformal coupling this can be done, thus achieving a second order fully debraided 
equation for a beyond-Horndeski model (see ref. [68] for an explicit debraided form in a more general case). 

We can think of the Jordan Frame version of Bekenstein’s theories as a simple playground to study the features of 
theories beyond Horndeski. In this sense they are analogous to the simple DBI Galileons discussed in section [ill B[ 
More general theories beyond Horndeski will not accept a simple Einstein frame formulation and will therefore not 
be as simple to unbraid. We expect that generalizations of Bekenstein theories will contain a richer mixing structure 
reflecting their more diverse phenomenology. Other theories beyond Horndeski have been proposed by Gleyzes, 
Langlois, Piazza and Vernizzi [31] . where it is also shown that no additional degrees of freedom are introduced and 
for which similar findings regarding kinetic mixing have been reported EDI321- In particular, the authors found terms 
describing the interactions between the scalar and derivatives of the matter energy density on perturbed cosmological 
backgrounds. See also refs. mm for further work on extensions beyond Horndeski. 


VI. CONCLUSIONS AND OUTLOOK 

In an era in which the alternatives to Einstein gravity are getting more and more complex, with an increasing 
level of mixing between the various degrees of freedom that build the theory, it is of fundamental importance to 
develop methods that allow the classification and facilitate the study of different models. This will not only lead to 
a better understanding of the physical content of a theory, but also can shed light on its properties, potential issues 
and observable consequences. 

In this work we have considered a method, covariant debraiding, to study the kinetic mixing between the scalar 
and tensor degrees of freedom in general, alternative theories of gravity. Our method consists in using contractions of 
the metric equations of motion to remove Ricci-curvature terms in the field equation of motion. This approach relies 
on the full equations of motion, and therefore allows to draw conclusions regardless of any approximation scheme, in 
a fully non-linear fashion with the extra advantage of being in the Jordan frame where the energy-momentum tensor 
of matter is covariantly conserved. Hence, the debraiding procedure provides a useful way to study the properties of 
ST theories and their interaction with matter as well as the stability of the scalar field rather directly. 

As an application of the method, we have extended covariant debraiding for the first time beyond the simplest 
examples and applied it to the study of quartic Horndeski theories in detail. These theories display a new set 
of mixing terms, which indicate new forms of coupling of the scalar field. The novel terms appearing in quartic 
Horndeski theories involve the contraction of second derivatives of the scalar with curvature which translate into 
contractions of and in the debraided equation. The procedure allows one to study how matter sources the 
scalar field despite both being coupled minimally, interacting only with the metric directly. 

General quartic theories also feature a coupling to the curvature that can not be removed by covariant debraiding. 
This term is given by the Weyl tensor, which is not algebraicly determined from the metric equations, contracted with 
the first and second derivatives of the scalar. It represents a novel form of interaction between the scalar field and 
space-time in the absence of matter. This form of coupling to the vacuum and the non-local nature of the Weyl tensor, 
determined by the global distribution of matter, might have relevant implications for Mach’s principle in ST theories, 
to be addressed in a follow up work. In addition, the debraided equation generally contain non-linear derivative terms 
without antisymmetric structure. These may introduce spurious solutions to the equations of motion, which can be 
nonetheless distinguished from the physical ones and disregarded. However, this problem is absent on sufficiently 
simple situations, such as cosmological backgrounds. 

Covariant debraiding singles out a particular subset of theories for which the debraiding equation neither contains 
the Weyl coupling, nor leads to spurious solutions. These theories generalize the quartic DBI Galileon, and coincide 
with the maximal set of Horndeski theories which accept an Einstein frame formulation (i.e. the kinetic term for the 
tensor degree of freedom can be written in the Einstein Hilbert form). This relation to Einstein gravity is behind 
the lack of Weyl coupling: as field redefinitions at the level of the action are equivalent to linear transformations of 
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the equations of motion, the occurrence of a term that can not be “rotated away” indicates the non-existence of a 
transformation to a non-mixed frame. 

Although the equations are far too involved for an analysis of general quartic theories, interesting conclusions can be 
easily drawn for the simple DBI-like models. Using the debraided equations we show that quartic DBI Galileons have 
gradient instabilities in the presence of matter with sufficiently large pressure. This would spoil the early universe 
predictions unless the energy scale that suppresses the coupling is very large or the theory is extended beyond the 
simplest case. In particular, a non-DBI quartic theory has a richer mixing structure, including terms which can 
stabilize the gradient instabilities. The debraided form of the equations that we present here can be used to design 
models with certain properties, by choosing the Horndeski functions to enhance a particular set of terms. 

Other features of quartic Horndeski theories are further clarified using their debraided formulation. One example 
is the screening of scalar forces by kinetic mixing. This mechanism was first investigated for DBI-like theories in the 
Einstein frame, and known as disformal screening mechanism. Our work implies that this effect is not exclussive of 
DBI-like/disformally coupled theories, but rather ubiquitous in quartic Horndeski theories. Moreover, more general 
theories might weaken the assumptions necessary for the screening mechanism to be effective. 

The mixing terms in the equations of motion are characteristic of each theory and carry the information about how 
the scalar degree of freedom interacts with matter. These terms grow in complexity in quartic theories (including the 
Weyl tensor coupling) and beyond, leading to non-linear mixing in quintic theories and derivative mixing in theories 
beyond Horndeski. The covariant debraiding procedure also provides a binary classification of models according to 
their kinetic mixing properties: Theories for which the covariant debraiding eliminates all instances of the curvature 
are stirred. This includes old-school theories, cubic and DBI-like quartic theories and some simple non-Horndeski 
theories. More generally, theories in which some residual curvature remains after covariant debraiding are shaken, 
including general quartic theories (featuring the Weyl tensor), quintic theories and general theories beyond Horndeski. 

There are other, potentially interesting avenues for further development. Our work has focused only on a very 
restricted form of debraiding, one in which locality and Lorentz invariance are manifest. One can use more general 
procedures to debraid the equations, e.g. by choosing a particular time slicing and solving for the highest time 
derivatives. This can be achieved (by explicit choice of coordinates or by an ADM decomposition), and is indeed 
necessary if one aims to numerically solve the general equations. As a step beyond, one may consider non-local 
debraiding by formally solving for Weyl tensor in terms of its propagation equation. These and other developments 
might shed light on a number of problems, such as the initial value formulation of modified gravity theories. 

The main lesson to be learned is that gravitational degrees of freedom become more fundamentally mixed the 
further away we go from Einstein’s theory of gravity. Any departure from the canonical kinetic term for the metric 
tensor necessarily introduces at least a scalar degree of freedom, which can then interact in a variety of ways, ranging 
from Brans-Dicke to beyond Horndeski. For both stirred to shaken theories, this increasing complexity reflects on 
the fundamental Lagrangian and gives important hints about its nature and dynamics. The study of kinetic mixing 
provides new ways to classify models and address their properties, valid across generations of ST theories, and will 
provide a useful tool in the understanding of alternative gravities. 
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Appendix A: Quartic Horndeski theories 


In this appendix we present the general equations for a Quartic Horndeski action. This kind of action assumes 
Gs{cj),X) = 0 = G§{<j>,X) while leaves the other two functions generic. We will first write the equations as derived 
from the variation of the Quartic Horndeski action @ and after present the detailed structure of the coefficients of 
the debraided equations that appears in (201. 
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1. Equations for quartic theory 

In this section we derive the both the metric and the scalar field equations in their explicit form. To sim¬ 
plify the expressions we will use the following notation for contractions [<f> n ] = (j>’ v , 4 = 

4>; fJ .a 1 4 l ’ ai - a2 •• (R) = 4)’ a 4>^R a p (if F a p has two indices), (W) = 4> ,a 4>’ l3 VV a ^4>’ 111 ' , (R(j>) = 4>’ a R a x4>’ XI3 4>^. 

The variation of the Lagrangians ([2]) and Q with respect to the scalar field, gives the following equation for the 
scalar field 


G 2 ,</> — 2 G 2 ^xX + (G 2 ,x — □</> — (G 2 ,xx + 2G 4)< ^x) 4 >,,J '4>’ v 4 >;ij, u 

+ 2 G^xx ( 2 (^ 2 ) - 2 {4>)[4>\ ~ {[4>] 2 - [0 2 ]) X) + G 4 ,xx ([</>] 3 - 3[</>][</> 2 ] + 2 [</> 3 ]) 

+ G^xxx (-2(<^ 3 ) + 2{cj> 2 )[(j)\ - (4>) ([0] 2 - [fi 2 ])) + 3G 4) 0x ([ <j )] 2 - [<^> 2 ]) 

+ G 4 ,^f? - 2 G^ X G^<!>,^ - 2G 4 ) 0x (2(i?) + RX) 

+ G 4>X x (-[</>]<!?) - \r[{4>) - [4>\X) + 2(W) - 2 (RtfiX + 2= 0 . 

(Al) 

The variation with respect to the metric gives 

GAG a p - l -T™ + - G 4 , 00 - 2 G^x[4>] - \ (kf - [</> 2 ]) - ^4,xi?) K4>V 

+ |(~G 2 + G 4 , 0 [</>] + ^G 4 , a ([</>] 2 - [cf 1 ] - (R)) - 2G iM X + \RG^xX 

—2 Gi t< f,x ((4 > ) + [<^]A) + G^xx {—(4 > )[4 > \ + (4> 2 ))) 9ap 

+ <Aa/3 (—G 4j 0 + G 4 > X.y(</>) — G 4 ,x[ 0] + 2G 4; 0xA r ) 

+ 2G 4j 0X {4^4 >r yp4 > a + '-ya4 > p) + Gk,X 4 > "t(S — RotpX + — Rp 1 4 )1 (j) a + ~R a ^4 >1 (j>p + Wa^p^cj)^ 

- G 4 ,XX (4 ,1 4 > p r] 4 ) 'ir,4a + 4> 1 4>a r] <i > 1 T ] 4>p - 4 )1 {4>'ip4 ) a + 4>'ya<f>p)[4>] + 4 }1 4 ) ia4 ) Pri4> V ) = 0 • (A2) 


2. Coefficients in the debraided equations 


In this appendix we report the explicit expressions for the coefficients of debraided equation (20 1 . These coefficients 
will be fixed once a particular model is chosen, i.e., a choice for the form of G 4 is made. 

The purely field dependent coefficients of the linear second order derivative are 


So 


Qt 

G(t) 


-(G 4 — 2G 4 j xA) 2 ^G 2 ,x( 3 G 4 + 8 G\ x X 2 — 4G 4 X(3G 4j x + G 4 i xxX)) + G 2 ( 4 G 4 i xX( 2 G 4i x 
+3G 4 ) xxA’) — G 4 (3G 4j x + 4G 4 i xxA)) — 3^4G 2 G 4 j< ^xX + G 4 (—3G 2 ^ 

+4G 4 j 0G 4 j 0xA + 4X(G 4 j xG 4 jfl j0 + 5G 2 ^ x X + 2 G 4 j a'xG 4 i 0 ^X — 4G 4 j xG 4 ^xA)) 
+ 8 G 4 j xA(G 2 ^ + A (—G 4 j a'G 4 j 0 ^ — 4G 2 0_y^ — 2G 4 j a'xG 4 0^A 


+ 2 G 4 | xG 4 | ^xA) 

(G 2 x + G 4 G 4 i xx)A 


(G 2 


3(G 4 - 2G 4 ,xA)(G 4 - G 4 ,xA) ’ 

G 4 G 4 i xx)(—3G 4 + 4G 4 j xA') 


4.X 


6 (G 4 - 2G 4 ,xA) 2 (G 4 - G 4 , x X) 


(A3) 

(A4) 

(A5) 
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3 G 2 ,aa (G 4 — 2G 4> xX) 2 — dGi^xiGA — 2G4 ; xA) 2 + ^G 4i xG 4>( j >( j > (—G4 + 2G 4: xX) 
+12 G 2 ^yA(—4G4 + 7 G 4j xA) + G 4 j x( — 3G2 ,xG 4 + G2G^\ — 3G 2 ^ + 4G2 i xG 4j xA) 

+G4 > 0Y ( — 24G4G 4j 0 + 36G4.yG4 j 0X) + G4 i A'a(G2G 4 + 12G4<^A( —G4 + 2G4 ; A'A) 

+4G2,xA(—2G4 + 3G 4> xA))) (G4 — 2G4 ; xA) 2 > 

c G‘i x + G 4 G 4j .ya 

T ~ 6 (G 4 - 2G 4 ,yA')(G 4 - G 4 ,a-A) ’ 

G\ x + G 4 G 4}X G 4tX x 
(T> = 6(G 4 - 2G 4 ,yA) 2 (G 4 - G 4 , a A) ■ 


(A6) 

(A7) 

(AS) 


The other kinetic couplings to matter/curvature read: 


Ct 

C(T) 

Cw 


G 4 ,x + G^xxX 
G 4 -G 4 , A-A ’ 

G\ x + G 4 G 4 ,aa- 

(G 4 -2G4 ,aA)(G 4-G 4 ,xA) ' 

2(Glx+G 4 G4,xx) 

G 4 - G 4 ,aA 


(A9) 

(A10) 

(All) 


Finally, there will be the generalizations of the (non-kinetic) conformal and disformal coupling, that we observed 
in KGB theories. 


Qt 

Q(t) 


Ga,4> + 2G4,a>A 
2(G 4 - 2G 4 ,aA) ’ 

4G 4 G 4j a> + G^x{G 4 ,4, — 6G 4j a<^A) 
2(G 4 - 2G 4 ,aA) 2 


The potential term is 


V = ^G2 1 (p{G 4 — 2G4 j aA)“ + G2 (—2G 4 G4 j( £ + G 4 j aA( 5G4 j 0 + 2G4 ^aA)) + A(— 2G2 t< j>x{G 4 
—2G 4 ^xX) 2 — &G 4> ^{G 4 — 2 G 4 ,aA)(G 4i 0 + 2 G 4 ^a'A) + G 2 ,a(G 4 G 4 ^ — 4G4 j a'G 4 j 0A 
-6G 4 G 4; 0 aA + 8 G 4 j aG 4j 0aA^ 2 ))^ (G 4 — 2G 4 . x X)~ 2 , 

The coefficients in front of the non-linear derivatives are 

Vb4 = ^G4(3G 4j 0a ^ 2G 4j 0aaA) + G4G4 j aA(G4 j aaA(— 13G 4i 0 + 3OG 4 j 0aA) 

- 4G4 i a( 3G4 i 0 - 11G4,0 aA + 4G4 i? jxaA 2 )) + G 2 4X X 2 (\2G 4} xxX{G 4 ^ - 2G 4 ^xX) 

+ G 4 i a( 11 G 4 ;( ^, — 30G4^aA + 8G 4t </,xxX 2 )) + G 4 (^G 4} xxX(3G 4t< i > — 10G 4 ^xX) 

+ G 4i x{3G 4t <i > — 21G 4 ^xX + 1 OG 4 j 0 xaA 2 ))^ (G 4 - 2G 4t xX)~ 2 (G 4 — G 4< xX )~ 1 , 


(A12) 

(A13) 


(A14) 


(A15) 


V 4D = (4G|G 4 ,<j>xx + G|(3G4 j aaG 4j 0 + 12G 4 ,aG 4 , 0 a + 22G 4t xxG 4 ^xX 

— 20G 4j aG 4 ^a'aA) + G 2 x A(48G 4j aaG4,0aA 2 + G 4) a(— SG^ + 3OG 4i 0xA 
— 16G 4j 0aaA 2 )) + G 4 G 4t x(—G 4i xxX(5G 4 ^ + 66G 4t< j,xX) + G 4 x(8G 4 ^ 

—38G 4 0aA + 32G4,^aaA 2 ))^ (G 4 — 2G4 i aA) _2 (G 4 - G 4j .yA) _1 , 


(A16) 
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V B5 = 


G 2 X + GaG^xX 


G 4 - Ga, x X ’ 

Vds = (3G|(G 4 ,^a - 2G^xxX) + GaGa, x X(SGaxxX{-Ga,4> + 12G 4 , 0 aA) 

+G 4j a(— 15G 4 ^ + 82G 4j( y\W — 48G 4i? iA.Y^ 2 )) + 4G 2 aA 2 (3G 4 i a'aA(G 4i 0 
— GG^xX) + G 4 j a(4G 4> 0 — ISGa^xX + QGa^xxX 2 )) + G 4 (—32G 4i axG 4i< ^a^ 2 
+3G 4 j a(G 4i ^ - UGa^xX + 1OG 4 ,0a.y^ 2 ))) (G 4 — 2G 4 i xA)“ 2 (G 4 — G 4 ) xAT ) _1 , 


(A17) 


(A18) 

The remaining terms are those we have labeled as “worrying terms”, as they might lead to spurious solutions on 
certain backgrounds (cf. Ill A). They read 


Wi = 


W 2 = 


(G 4i a + G 4 G 4 i aa-)(G4,x + 2G 4) xxX) 
3(G 4 -2G 4 iX A)(G 4 -G 4 , a -X) 

(G 2 x + G 4 G 4j aa ) 2 


W D2 = - 


3(G 4 - 2G 4 ,. y X) 2 (G 4 - G 4 , x X) ’ 

[G\ x + G 4 G 4i aa)(4G 4 G 4j 0.y + G 4i a(G 4j ^ — GGa^xX)) 

3(G 4 -2G 4 ,aA') 2 (G 4 -G 4i aX) 


(A19) 

(A20) 

(A21) 


a. DBI-like quartic theories 


We report here the detailed expressions for the DBI-like quartic theories introduced in EH 

Go = G 2 , A + 

8G2 A 4 (A(</>)) 3 VX 2 + M|pY 8 :D 4 A(-13A(A'(</>)) 2 + 4A 4 X> 2 A"(</>) 
+4(A(</>)) 2 (-5A/| 1 X 4 (A , (<^)) 2 + 2A 4 X> 2 A (G 2 A 4 X> + 2 M^X 2 A' 

+2A 4 A(<£)X> 2 14M|jA 3 (A'(</>)) 2 + A 4 !? 2 (G 2 A 4 £> + SM^A 2 A' 

x (2A 8 P 3 (A 4 X> 2 + 2A(<f>)X) 


So — ~ G 2 ,aa + 

4(A(</>)) 3 (2G 2 ,xA 8 £> 3 A 2 - 3Ml, l X 4 (A'(4 >)) 2 ) + 2M^A 12 X> 6 (-10A(A'(</>)) S 
+A 4 X> 2 A"(^)) + 4A 4 (A(^)) 2 X> 2 a(-13A/| 1 A 2 (A , (</>)) 2 + 2 A 4 X> 2 (G 2 ,a A 4 !? 
+M 2 j A A" (</>))) + A 8 A(^)T> 4 (-59M 2 1 A 2 (A , ( < />)) 2 + 2A 4 V 2 (G 2 , x A 4 V 
+4M| 1 AA//(0))) (2A 12 X> 5 (A 4 X> 2 + 2A(<j))Xy 


V = G 2 ,*-2G ; 




X 


XA'((f)) ^A 4 A(0)'D 2 a(-24M 2 1 A 3 (A , (0)) 2 + A 4 D 2 (A 4 D(7G 2 + 10G 2 ,aA) , 
-24A/| 1 A 2 A"((^))) + A 8 T> 4 (-9M| 1 A 3 (A , (^)) 2 + A 4 V 2 (A 4 V{2G 2 + 5G 2 ,aA) 

- 2(A(^)) 2 A 2 (-6A// 2 1 A 3 (A , (0)) 2 + A 4 V 2 (A 4 V(G 2 + 4G 2 ,aA) 


—9Afp 1 A 2 A" 


—6Afp 1 A 2 A" 


(a 12 X> 5 (A 4 T> 2 + 2A(</))A) : 


(A22) 


(A23) 


(A24) 
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where we have defined 


V4 d = 


2 M&AMA'ty) 
A 8 T > 3 


V4 B = - 


M|j(3A 4 D 2 - 2 A{<j>)X)A'(<l>) 


2A S V 3 


X (3A 4 !? 2 + 2A{<j))X)A'((t>) 


Q(t) - 


Qt ~ 2A 8 T> 4 + 4A 4 A(</>)2? 2 A' ’ 

(4A 8 T > 4 + 9A 4 A(cj))'D 2 X + 6(A(0)) 2 A 2 ) A'(</>) 
2A 4 (A 4 V 3 + 2 A{cj))VX) 2 

m 


Ct — 


A 4 V 2 


V = (1- 2A(</>) X/A 4 ) 1/2 


(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 


Appendix B: Existence of an Einstein frame 


In section [TV A| we have discussed how, the presence of the Weyl tensor in the equation of motion for the scalar field 
can be related to the lack of a held transformation able to cast the kinetic term for the metric into its standard Einstein- 
Hilbert form. In this appendix we provide additional arguments in support of this statement. More specifically we 
examine the requirements for a Held redefinition to be able to cast a quartic theory as Einstein-Hilbert plus a coupling 
to matter. In what follows we enumerate different possibilities, and argue that they are implausible to produce the 
desired outcome. 

1. If we use a scalar Held redefinition only we need G^X, <j>) — > M 2 /2 = constant. However, we can see already in 

the old-school case that demanding that G 4 (</>(£)) = constant implies = 0 (by taking the derivative wrt 

the field). This requires that either G 4 is constant or that the Jacobian of the field transformation is degenerate. 
It is hard to imagine how adding X dependence would help with this. 

2. If we use a metric redefinition we note that 


(a) The Riemann tensor transforms as 


R a p^ = R a ^ + w [ll ic a v]p + 




(Bl) 


where /C Q w = f“ M „ — T a [30]. So any factor that cancels the coefficient of R has to come from either 
\f-~g or from . 

(b) Using X dependent disformal transformations can neutralize all the dependences of Gi(X : (j>), but non- 
Horndeski terms are generated in the transformation |46j . These do not lead to an unhealthy theory (see 
appendix [C]), but does not provide a canonical kinetic term for gravity either. Restricting to ^-dependent 
disformal transformations prevents the non-Horndeski terms from appearing, but then only DBI-like forms 
of G 4 can be canceled. 


(c) If we try to use higher derivatives for the metric redefinition ■ ■ ■ ...) then we will introduce these 

higher derivatives in G 4 unless every dangerous term cancels. This can be seen from the structure of (Bl I: 
The derivative dependence from g a P and \J—g will appear in G 4 , unless the connection terms produce 
a coefficient proportional to R that cancels all the field dependedences in G 4 , with the exception of a 
constant term and without introducing other curvature terms. Moreover, one would in general introduce 
higher derivatives via the second term in [Bl[ that also need to cancel. 


The above reasoning just argues the implausibility of finding a local field redefinition that transforms the kinetic term 
for quartic theories in a canonical form. It does not constitute a mathematical proof, but complements the discussion 
presented in IV A In particular, our discussion does not apply to non-local field redefinitions. 
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Appendix C: Beyond Horndeski: Examples 

In this appendix we briefly present the simplest theories beyond Horndeski, the pure conformal and pure disformal 
theories with derivative dependences. These theories can be formulated as Einstein Hilbert plus a coupling to matter 
via a field redefinition. In this sense they are similar to the DBI-Galileons, and they provide toy examples of kinetic 
mixing in theories beyond Horndeski. 


1. Pure Conformal Theory 

In the relatively simple case of an X-dependent conformally coupled theory in the Jordan frame 

V~9 


C ° 16t tG 

one can write the equations of motion as 


(H 2 R + + yj—g + C m ) , 


(Cl) 


^T K ) + §±T, k-i^ = o. 

where the kinetic mixing factor for a conformal relation reads 

_ 8 nGQ,xT 

K - ci - 2 n >x x' 


(C2) 


(C3) 


with T = g^ v (T$ v + TJH). The derivative mixing enter ( in the Jordan frame) as a way to remove the higher order 
derivatives of the scalar field in this type of theories. 


2. Pure Disformal Theory 

Let us now consider the transformation of the Einstein-Hilbert Lagrangian under a disformal transformation de¬ 
pending on field derivatives 


g^u — g/iv T , (C4) 

whose associated barred connection is characterized by the following lC a ^ v = Yf lv — T“ y tensor \SU 

IC a - (log H),*f ~ 7V“<*>] , (C5) 

We will denote 7 2 = U 7 2 = B/( 1 - 2 BX), (<&) = , (4> 2 ) = [$] = □<£, [•L 2 ] = and 

(R) = (j)’^. The first term was already present in the only field-dependent disformal transformation, while the 
last two arise from the derivatives of B. The bulk contribution to the Einstein-Hilbert action can be computed using 
Eq. (38) of gD], yielding 

= V=g(^R- R'y(R) - b(b + h,xX) 7 3 (($)[$] - ($ 2 ))^) . (C6) 


This simple result is due to the particular tensor structure of the K, tensor, which makes that only the first term 
in Eq. (C5) contributes to (all other contractions are proportional to a contraction of <j> ,a <j>’^ with a 

tensor antisymmetric on the indices ag). The addition of a surface term V M (/(X)(0 ,#i [<j>] — </> ;/i “</> a )) to (C 6 ) gives 


V=$R = V=g(^R-{'yB + f){R) + fm 2 -[* 2 ]) 

- (B(B + b, a-x ) 7 3 + f. x ) (($)[$] - ($ 2 ))). 


(C7) 


It is not possible to write down the above action in the canonical Horndeski form, as that would require to simulta¬ 
neously satisfy 
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1 . / = (I/ 7 ), x = —7 B — 7 BxX, to have the right coefficient of [$] 2 — [$ 2 ], 

2. f\x = —B(B + B^xX'j'y 3 , to kill the last term and 

3. / = —7 B, to kill the (R) term. 

The last choice gives the simplest form for the action 


S==gR = y /^(±R- Bl m 2 - [<h 2 ]) - 7 Sx «*>[*] - (<& 2 ») 


(C 8 ) 


Similarly to the conformal theory, the equations of motion for the theory (C 81 contain an implicit constraint that 


allows one to remove the higher derivatives in the equations of motion. In this case the kinetic mixing factor reads 


r D 


K 


8'rrG'y 


1 + 2B yX X 2 


(C9) 


and also enters the field equation as X a (T x </>’“) [30] . 
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